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PREFACE 


The  numerical  model  investigation  of  the  San  Francisco  River  through 
Clifton,  AZ,  reported  herein,  was  conducted  at  the  US  Army  Engineer  Waterways 
Experiment  Station  (WES),  at  the  request  of  the  US  Army  Engineer  District, 

Los  Angeles  (SPL). 

This  investigation  was  conducted  during  the  period  May  to  December  1988 
in  the  Hydraulics  Laboratory  of  WES  under  the  direction  of  Mr.  Frank  A. 
Herrmann,  Jr.,  Chief  of  the  Hydraulics  Laboratory,  Mr.  Marden  B.  Boyd,  Chief 
of  the  Waterways  Division  (WD) ,  and  Mr.  Michael  J.  Trawle,  Leader  of  the  Math 
Modeling  Group  (MMG) .  Mr.  William  A.  Thomas,  WD,  provided  general  guidance 
and  review.  The  project  engineer  was  Mr.  Ronald  R.  Copeland,  MMG.  Authors  of 
this  report  were  Mr.  Copeland  and  Mr.  Thomas.  Technical  assistance  was  pro¬ 
vided  by  Ms.  Brenda  L.  Martin,  MMG. 

During  the  course  of  this  study,  close  working  contact  was  maintained 
with  Mr.  Rene  Vermeer en  and  Ms.  Rachel  Korkos  of  the  Hydraulics  Section,  SPL, 
who  served  as  coordinating  engineers,  providing  required  data  and  technical 
assistance. 

Acting  Commander  and  Director  of  WES  during  preparation  of  this  report 
was  LTC  Jack  R.  Stephens,  EN.  Technical  Director  was  Dr.  Robert  W.  Whalin. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement 
(metric)  units  as  follows: 

_ Multiply _ 

cubic  feet 
cubic  yards 
feet 
inches 

miles  (US  statute) 
square  miles  (US  statute) 
tons  (2,000  pounds,  mass) 


used  in  this  report  can  be 

_ _ 

0.02831685 

0.7645549 

0.3048 

2.54 

1.609347 

2.589998 

907.1847 


converted  to  SI 

To  Obtain 
cubic  metres 
cubic  metres 
metres 
centimetres 
kilometres 
square  kilometres 
kilograms 
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EVALUATION  OF  BANK  EROSION  POTENTIAL,  CLIFTON,  ARIZONA 


Numerical  Model  Investigation 


PART  I :  INTRODUCTION 

Description  of  Study  Area 

1.  The  town  of  Clifton  is  located  in  mountainous  terrain  in  southeast¬ 
ern  Arizona  (Figure  1).  The  San  Francisco  River  flows  through  the  town  toward 
its  confluence  with  the  Gila  River.  The  river's  drainage  area  above  Clifton 
is  approximately  2,800  square  miles.*  Elevations  in  the  basin  range  from 
9,200  ft**  near  the  headwaters  to  3,500  ft  near  Clifton.  A  single  levee  is 
proposed  to  protect  south  Clifton  for  a  one-half-mile  reach  adjacent  to  the 
San  Francisco  River  beginning  immediately  below  the  Southern  Pacific  Railroad 
and  Highway  666  Bridges  (sta  568+00)  and  extending  downstream  to  Ward  Canyon 
at  sta  534+00  (Figure  2) .  The  project  reach  extends  upstream  through  north 
Clifton  for  a  total  project  length  of  about  2.3  miles.  The  project  upstream 
of  the  railroad  bridge  would  consist  of  housing  relocation  and  floodproofing 
but  no  improvements  to  the  existing  channel.  The  project  is  designed  for  the 
125-year-frequency  flood  of  145,000  cfs.  The  Standard  Project  Flood  (SPF)  is 
167,000  cfs. 

2.  Existing  channel  improvements  through  Clifton  include  masonry  flood- 
walls  and  dumped  molten  slag  deposits  from  local  copper  smelters  (Figure  2). 
These  bank  protection  measures  were  constructed  between  1906  and  1920.  The 
right  descending  bank  downstream  from  the  railroad  and  highway  bridges  has 
slag  deposit  bank  protection  for  1,200  ft  and  then  a  nearly  vertical  bluff 
wall  composed  of  rock  or  slag  through  the  rest  of  the  study  reach.  Downstream 
from  the  bridges,  on  the  left  descending  bank,  a  masonry  wall  extends  for 
about  1,100  ft  until  it  intersects  a  slag  deposit  of  uncertain  continuity.  An 
investigative  study  (US  Army  Engineer  District  (USAED) ,  Los  Angeles,  1988) 
concluded  that  this  slag  deposit  was  probably  continuous  between  river 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  found  on  page  3. 

**  All  elevations  (el)  cited  herein  are  in  feet  referred  to  the  National  Geo¬ 
detic  Vertical  Datum  (NGVD). 
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ARIZONA 


Clifton,  AZ,  location  map 


Clifton,  AZ,  project  features 


sta  535+00  and  sta  562+00  except  for  a  150-ft-long  section  near  sta  541+00. 
However,  the  existence  of  the  slag  deposit  between  sta  548+00  and  sta  562+00 
is  uncertain  because  the  area  is  covered  by  fill  material.  The  existence  of 
the  slag  deposit  was  inferred  based  on  conversations  with  local  residents  and 
dark  lineaments  on  aerial  photos.  The  deposit  was  reportedly  dumped  from  nar¬ 
row-gaged  tracks  laid  along  the  riverside  edge  of  mine  tailing  piles.  The 


condition  and  configuration  of  the  deposit  varies.  Thickness  varies  between  3 
and  25  ft;  height  between  5  and  18  ft;  toe  elevations  from  below  the  waterline 
to  5  ft  above  the  waterline;  slopes  from  near-vertical  to  IV  on  4H;  and  foun¬ 
dation  material  from  loose  silt  to  hard  conglomerate.  The  slag  was  found  to 
be  generally  competent,  with  zones  of  deterioration  and  slumping. 

3.  The  alignment  of  the  proposed  levee  is  set  back  about  100  ft  from 
the  existing  main  river  channel.  The  100-ft  setback  was  established  based  on 
hydraulic  calculations  that  predicted  low  velocities  adjacent  to  the  levee 
toe,  on  historical  flood  indications  that  the  riverbank  is  stable,  and  on  the 
geotechnical  analysis  of  the  bank  soils.  The  levee  alignment  was  also  coordi¬ 
nated  with  local  interests.  The  proposed  earth  levee  includes  grouted  stone 
side-slope  protection  extending  5  ft  below  existing  ground  elevation.  The 
levee  is  2,500  ft  long  with  IV  on  2.5H  side  slopes.  The  levee  ranges  in  height 
from  7  to  17  ft,  with  a  minimum  3  ft  of  freeboard.  Two  manually  operated 
floodgates  are  proposed  at  the  intersection  of  the  levee  alignment  and  the 
railroad  and  highway.  These  structures  include  swing  and  roller  steel  flood¬ 
gates  15  ft  high,  and  concrete  retaining  walls  having  a  maximum  height  of 

18  ft. 

Purpose  and  Approach 

4.  The  purpose  of  this  study  was  to  evaluate  the  potential  for  bank 
erosion  adjacent  to  the  proposed  levee.  This  was  accomplished  indirectly  by 
identifying  potential  for  degradation  and  meandering  in  the  main  river  chan¬ 
nel.  Degradation  potential  was  evaluated  using  the  TABS-1  one-dimensional 
numerical  model.  Meandering  potential  was  evaluated  by  reviewing  historical 
river  data  and  by  application  of  generalized  empirical  regression  equations 
for  channel  width. 

5.  The  numerical  model  study  evaluated  the  potential  magnitude  of  scour 
and  deposition  in  the  study  reach  for  the  SPF,  and  for  a  flood  with  the  same 


peak  discharge  as  the  flood  of  record.  Model  adjustment  was  based  on  the  as¬ 
sumption  chat  the  existing  channel  is  essentially  stable.  The  effect  of  back- 
wat^  upstream  from  the  railroad  and  highway  bridges  on  sediment  transport  and 
bed  change  was  determined. 

6.  Historical  river  data,  including  river  profiles,  aerial  photos,  and 
topographic  maps,  were  reviewed.  The  historical  data  were  supplemented  by  a 
field  reconnaissance  to  observe  existing  channel  conditions.  These  data  were 
used  to  evaluate  the  historical  stability  of  the  river. 

7.  Generalized  empirical  regression  equations  that  predict  stable  chan¬ 
nel  width  were  used  to  evaluate  the  potential  for  bank  erosion.  Due  to  the 
generalized  nature  of  these  equations,  their  regional  applicability,  and  the 
uncertainty  related  to  the  influence  of  bank  stability,  results  from  this  type 
of  analysis  are  valuable  only  in  the  sense  that  they  provide  knowledge  of 
probable  trends  in  channel  geometry  change. 
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PART  II:  INPUT  DATA  DESCRIPTION 

Description  of  Numerical  Model 


8.  The  TABS-1  one-dimensional  sedimentation  program  was  used  to  develop 
the  numerical  model  for  this  study.  Development  of  this  computer  program  was 
initiated  by  Mr.  William  Thomas  at  the  US  Army  Engineer  District,  Little  Rock, 
in  1967.  Further  development  at  the  Hydrologic  Engineering  Center  (HEC)  by 
Mr.  Thomas  produced  the  widely  used  HEC-6  generalized  computer  program  for 
calculating  scour  and  deposition  in  rivers  and  reservoirs  (HEC  1977).  Addi¬ 
tional  modification  and  enhancement  to  the  basic  program  by  Mr.  Thomas  at  the 
US  Army  Engineer  Waterways  Experiment  Station  (WES)  led  to  the  TABS-1  program 
currently  in  use.  TABS-1  is  considered  experimental  in  that  it  is  not  docu¬ 
mented  to  the  point  that  it  can  be  made  available  for  general  use,  but  can  be 
made  available  by  special  request.  The  program  produces  a  one-dimensional 
model  that  simulates  a  series  of  steady-state  discharge  events  and  their 
effect  on  the  sediment  transport  capacity  at  cross  sections  and  the  resulting 
degradation  or  aggradation. 

Channel  Geometry 


9.  Channel  geometry  for  the  numerical  model  was  developed  from  a  HEC-2 
backwater  model  prepared  by  the  Los  Angeles  District  in  December  1985,  and 
adjusted  in  February  1986.  This  model  included  the  proposed  levee  and  other 
design  features  and  extended  from  sta  514+42  upstream  to  sta  653+20.  Cross- 
section  elevations  in  the  channel  were  taken  from  a  1977  survey  of  the  chan¬ 
nel,  and  overbank  elevations  from  a  1939  topographic  map.  The  TABS-1  sediment 
model  does  not  have  a  bridge  analysis  capability.  Therefore,  backwater 
effects  from  the  highway  and  railroad  bridges  were  accounted  for  by  setting 
the  water-surface  elevation  upstream  from  the  bridges  based  on  a  rating  curve 
developed  from  the  HEC-2  backwater  calculations.  The  upstream  boundary  of  the 
numerical  model  at  sta  653+20  was  artificially  extended  5,000  ft  to  remove  the 
boundary  from  backwater  effects  caused  by  the  constricted  channel  sections  in 
the  vicinity  of  the  Patterson  Addition  Access  Bridge  near  sta  626+50  (Fig¬ 
ure  1).  This  was  accomplished  in  the  model  by  repeating  the  upstream  cross 
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section  five  times  at  intervals  of  1,000  ft  and  increasing  bed  elevations 
based  on  a  slope  of  0.0033. 


Hydrographs 

10.  The  discharge  hydrographs  were  simulated  in  the  numerical  model  by 
a  series  of  steady-state  events.  The  duration  of  each  event  was  chosen  such 
that  changes  in  bed  elevation,  due  to  deposition  or  scour,  did  not  signifi¬ 
cantly  change  the  hydraulic  parameters  during  that  event.  At  relatively  high 
discharges,  durations  need  to  be  short.  At  low  discharges,  the  time  interval 
may  be  extended.  Time  intervals  between  15  minutes  and  2  hours  were  used  in 
this  study. 

11.  A  hydrograph  simulated  by  a  series  of  steady-state  events  of  vary¬ 
ing  durations  is  called  a  histograph.  The  SPF  histograph  used  in  the  numerical 
model  investigation  was  based  on  a  hydrograph  developed  by  the  Los  Angeles 
District  (USAED,  Los  Angeles,  1986).  This  flood  had  a  peak  discharge  of 
167,000  cfs.  The  histograph  is  shown  in  Plate  1. 

12.  A  flood  with  a  peak  discharge  equal  to  that  of  the  flood  of  record 
(90,900  cfs  in  1983)  was  also  tested.  This  histograph  was  assumed  to  have  the 
same  duration  as  the  SPF,  and  each  discharge  event  was  reduced  by  the  ratio  of 
the  flood  peaks  (0.544).  This  histograph  is  shown  in  Plate  2. 

Downstream  Water-surface  Elevation 


13.  The  water-surface  elevation  rating  curve  at  the  downstream  boundary 
(sta  514+42)  was  based  on  the  HEC-2  calculations  by  the  Los  Angeles  District 
(Plate  3).  The  rating  curve  at  the  downstream  boundary  was  based  on  the 
slope-area  option  in  the  HEC-2  code.  A  rating  curve  upstream  from  the  rail¬ 
road  and  highway  bridges  (sta  566+80)  was  based  on  District  desk-top  calcula¬ 
tions  using  orifice  and  weir  equations  at  high  flow  and  the  HEC-2  special 
bridge  routine  at  low  flows  (Plate  4). 

Bed  Material  and  Overbank  Gradations 


14,  Size-class  characteristics  of  the  bed  material  and  the  overbank 
were  developed  from  1985  test  trench  data  (USAED,  Los  Angeles,  1986).  These 
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data  were  collected  as  part  of  a  geologic  and  soils  investigation  to  determine 
suitable  borrow  sites,  foundation  conditions,  and  soil  design  values  for  the 
levee  project.  Locations  of  the  test  trenches  are  shown  in  Plate  5. 

15.  Ten  test  trenches  were  dug  in  the  riverbed  to  depths  of  12  ft. 
Gradations  of  samples  from  these  test  trenches  were  determined  by  sieve  analy¬ 
ses.  Although  boulders  with  diameters  as  large  as  2  ft  were  observed  in  the 
trenches,  all  material  larger  than  3  in.  was  excluded  from  the  samples.  No 
field  estimate  of  the  percentage  of  this  coarse  material  in  the  bed  was  made. 
In  the  numerical  model,  8  percent  of  the  bed  material  was  arbitrarily  assumed 
to  be  greater  than  3  in.  An  analysis  of  the  test  trench  data  revealed  that 
all  but  two  of  the  17  sample  gradations  could  be  identified  as  either  sandy 
silt  or  gravelly  sand.  Interestingly,  samples  from  a  given  site,  but  at  dif¬ 
ferent  depths,  could  fall  into  both  groups.  This  vertical  variation  in  bed 
material  indicates  active  bed  movement  either  due  to  bedforms  or  shifting 
gravel  bars  and  sinuosity  patterns.  Average  gradation  curves  for  the  two 
types  of  samples  and  their  ranges  are  shown  in  Plate  6.  Also  shown  are  the 
two  intermediate  samples  that  did  not  fall  into  either  category.  The  average 
gradation  for  the  gravelly  sand  material  is  considered  most  representative  of 
the  active  layer  of  bed  material  during  major  flow  events  and  was  used  in  the 
numerical  model.  A  sensitivity  study  was  conducted  to  determine  the  effect  of 
using  the  coarser  and  finer  limits  of  the  gravelly-sand  gradations. 

16.  Lateral  and  longitudinal  variations  in  the  bed  material  gradation 
can  significantly  affect  scour  and  deposition  in  gravel-bed  rivers.  These 
variations  can  be  especially  important  at  channel  constrictions,  expansions, 
and  bends.  However,  the  San  Francisco  River  data  are  insufficient  to  define 
either  the  lateral  or  longitudinal  variation  in  bed  material  gradation. 

17.  On  the  overbank,  five  trenches  were  dug  along  the  proposed  levee 
alignment  to  depths  of  12  ft.  These  data  indicate  that  the  overbank  is  com¬ 
posed  primarily  of  mine  tailings  consisting  of  stratified  layers  of  clay, 
clayey  silt,  silty  clay,  silt,  silty  sand,  sandy  gravel,  and  gravelly  sand. 
Detailed  size-class  analysis  was  not  conducted,  but  general  size  descriptions 
show  that  the  silt-clay  percentage  in  the  overbank  varies  between  18  and 

100  percent.  The  sand  and  gravel  layers  give  the  overbank  a  high  potential 
for  undermining  and  bank  failure.  The  upstream  500  ft  of  the  proposed  levee 
would  rest  on  a  fairly  erosion-resistant  conglomerate  formation. 


Roughness 

18.  Channel  roughness  In  the  TABS-1  numerical  model  was  based  on  re¬ 
sults  from  the  HEC-2  backwater  study.  Manning's  roughness  coefficients  in  the 
backwater  model  had  been  adjusted  to  reproduce  high-water  stages  from  histori¬ 
cal  flood  events.  Channel  roughness  coefficients  varied  between  0.025  and 
0.030.  The  channel  roughness  coefficient  at  the  upstream  cross-section 
(653+20)  in  the  HEC-2  model  was  changed  from  0.050  to  0.025  in  the  TABS-1 
model  to  be  consistent  with  downstream  values.  Overbank  roughness  coeffi¬ 
cients  varied  between  0.045  and  0.060.  Major  obstructions  in  the  overbank 
were  accounted  for  by  inclusion  in  the  cross-section  geometry. 

Sediment  Inflow 

19.  Between  1982  and  1984,  suspended  sediment  samples  were  collected  at 
the  Clifton  gage,  which  is  located  about  one  mile  downstream  from  the  town  at 
the  sewage  treatment  plant  (sta  515+00).  The  maximum  discharge  at  which  sam¬ 
ples  were  obtained  was  2,000  cfs.  Particle  size  analyses  were  not  performed 
on  these  samples  so  silt  and  clay  percentages  are  unknown.  The  data  are  use¬ 
ful  only  for  general  comparisons  with  calculated  bed-material  loads. 

20.  Sediment  inflow  for  the  TABS-1  numerical  model  was  estimated  by 
assuming  equilibrium  transport  potential  at  the  upstream  boundary.  It  was 
necessary  to  extend  the  model's  upstream  boundary  5,000  ft  due  to  backwater 
created  by  the  constriction  near  the  Patterson  Addition  Access  Road  Bridge 
(sta  626+50).  Due  to  the  uncertainty  related  to  the  sediment  inflow,  sensi¬ 
tivity  calculations  were  made  with  sediment  inflows  greater  and  less  than 
estimated  average  values. 

Transport  Function 

21.  Four  transport  functions  were  evaluated  for  use  in  the  numerical 
model.  These  were:  Meyer-Peter  and  Muller  (1948)  bedload  formula,  Laursen- 
Madden  (HEC  1977)  bed  material  load  formula,  a  combination  of  the  Schoklitsch 
(Schulits  1935)  bedload  formula  and  the  Toffaleti  (1966)  bed-material  load 
formula,  and  a  combination  of  the  Meyer-Peter  and  Muller  and  Toffaleti  formu¬ 
las.  Hydraulic  parameters  for  the  evaluation  were  taken  from  fixed-bed 
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backwater  runs  at  the  upstream  cross  section.  Results  of  calculations  using 
these  functions  are  compared  to  measured  suspended  data  in  Plate  7.  These 
data  are  comparable  only  in  a  general  sense  because  the  measured  data  include 
silts  and  clays  and  excludes  bedload,  while  the  calculated  data  include  total 
sand  and  gravel  transport  but  exclude  silts  and  clays.  The  combination  of 
Toffaleti  and  Meyer-Peter  and  Muller  formulas  was  selected  because  it  appeared 
to  give  reasonable  results  for  sand  transport  and  was  capable  of  transporting 
the  larger  gravel  sizes.  Sediment  inflow  rating  curves  for  each  size  class 
used  in  the  numerical  study  are  shown  in  Plates  8  and  9. 
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PART  III:  STUDY  RESULTS 


Average-Annual  Flood 

22.  To  get  a  general  idea  of  aggradation  and  degradation  potential  in 
the  study  reach,  the  numerical  model  was  used  to  simulate  a  steady-state  dis¬ 
charge  of  7,500  cfs  (2.33-year  frequency)  for  a  duration  of  five  days.  This 
simulation  indicated  a  degradation  trend  downstream  from  the  railroad  and 
highway  bridges.  The  model  also  indicated  that  the  riverbed  through  Clifton 
is  highly  dynamic,  with  sediment  transport  potential  varying  significantly 
with  distance  along  the  channel.  Due  to  the  dynamic  characteristics  of  the 
sediment  movement,  it  was  decided  that  the  SPF  hydrograph  should  be  simulated 
to  determine  the  effect  of  duration  on  scour  and  deposition  in  the  study 
reach. 


Standard  Project  Flood 

23.  The  SPF  histograph  was  run  with  1977  channel  geometry  and  equilib¬ 
rium  sediment  inflow.  Insufficient  data  were  available  to  adjust  or  circum¬ 
stantiate  the  numerical  model,  therefore  results  can  be  used  only  in  a  quali¬ 
tative  sense.  This  qualitative  evaluation  is  valuable  in  determining  the 
potential  for  bank  erosion  adjacent  to  the  proposed  levee  alignment.  Sections 
that  show  significant  scour  or  aggradation  are  indications  of  a  section  with 
high  erosion  potential. 

24.  The  model  showed  that  the  riverbed  through  Clifton  is  highly  dy¬ 
namic  during  the  passage  of  the  SPF.  Thalweg  profiles  at  the  peak  and  at  the 
end  of  the  simulation  are  compared  to  the  initial  thalweg  in  Figure  3.  This 
figure  shows  scour  on  the  order  of  5  ft  downstream  from  the  railroad  and  high¬ 
way  bridges  at  the  peak  of  the  flood.  This  scour  trend  does  not  continue 
throughout  the  leveed  reach  because  of  the  constriction  at  sta  521+80.  This 
constriction  creates  a  backwater  that  prevents  degradation  at  discharges 
greater  than  40,000  cfs.  The  backwater  effect  disappears  at  lower  discharges 
which  helps  explain  the  incised  channel  observed  in  the  field.  The  constric¬ 
tion  itself  is  subject  to  significant  degradation  which  results  in  deposition 
downstream.  This  condition  does  not  make  for  a  favorable  downstream  boundary 
in  the  model.  Downstream  control  should  be  investigated  in  detail  to 
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determine  if  there  is  invert  control  at  the  constriction.  Bed  changes  during 
the  course  of  the  SPF  simulation  at  five  cross  sections  in  the  leveed  reach 
are  shown  in  Plates  10-14. 

25.  The  model  results  upstream  from  the  leveed  reach  are  also  interest¬ 
ing.  The  San  Francisco  River  is  constricted  through  the  town  of  Clifton,  and 
as  a  result,  significant  scour  occurs  at  the  first  sections  inside  the  im¬ 
proved  reach  of  river.  The  model  showed  up  to  13  ft  of  scour  near  the  SPF 
peak.  Unless  the  Patterson  Addition  Access  road  bridge  is  built  on  bedrock  or 
has  extensive  pile  footings,  its  integrity  is  questionable  during  a  major 
event.  In  general,  the  entire  reach  of  river  through  Clifton  is  subject  to 
degradation  during  a  major  flood.  This  is  due  to  the  increased  velocities 
caused  by  the  constriction  through  town  and  the  backwater  created  by  the 
constriction  upstream  from  town  that  reduces  sediment  transport  potential  to 
the  improved  reaches.  Backwater  from  bridges  and  local  constrictions  serve  to 
create  local  areas  of  deposition  and  scour. 

Sensitivity  Analysis 

26.  Sensitivity  of  the  model  to  sediment  inflow  was  tested  by 
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increasing  and  decreasing  the  sediment  inflow  for  each  size  class  by  a  con¬ 
stant  ratio  of  2.0  and  0.5,  respectively.  Comparisons  of  the  total  sand  and 
gravel  inflow  curves  to  the  measured  data  are  shown  in  Figure  4.  Calculated 
thalweg  elevations  at  the  peak  of  the  SPF  with  the  different  sediment  inflow 
assumptions  are  compared  in  Figure  5.  Calculated  thalweg  elevations  at  the 
end  of  the  SPF  are  compared  in  Figure  6.  These  comparisons  show  that  bed  ele¬ 
vations  in  the  leveed  reach  are  not  significantly  affected  by  this  range  of 
sediment  inflow  assumptions. 

27.  Additional  sensitivity  studies  were  conducted  to  determine  the 
effect  of  the  bed  material  gradation  on  degradation  and  aggradation  in  the 
study  reach.  The  average  bed  material  gradation  from  the  test  trench  data  had 
a  D,_g  =  3.5mm.  A  coarser  gradation,  =  9mm,  and  a  finer  gradation,  = 
1.5mm,  were  obtained  from  the  limits  of  the  gradation  data  for  gravelly-sands 
from  the  test  trenches  (Plate  6).  Results  of  the  sensitivity  test  are  shown 
in  Figure  7.  The  bed  material  gradation  did  not  significantly  affect  scour  or 
deposition  in  the  study  reach.  The  sensitivity  test  does  not  address  the  pos¬ 
sible  effects  of  longitudinal  variation  in  the  bed  material  gradation. 

28.  A  sensitivity  test  was  conducted  to  determine  the  calculated  scour 
and  deposition  for  a  flood  event  similar  to  the  1983  flood  of  record.  Results 
of  the  simulation  are  shown  in  Figure  8,  and  compared  to  SPF  results  in 
Plates  15  and  16.  The  magnitude  of  scour  for  the  1983  simulated  flood  was 
similar  to  the  SPF  scour  in  the  leveed  reach. 

Historical  Analysis 

29.  Historical  accounts  of  numerous  floods  show  that  the  alignment  of 
the  main  river  channel  has  remained  relatively  straight,  and  the  banks  essen¬ 
tially  stable,  through  the  proposed  levee  reach.  Bank  erosion  along  the  left 
descending  bank  adjacent  to  the  proposed  levee  was  reported  during  interviews 
with  local  residents  after  the  1983  flood  (USAED,  Los  Angeles,  1986).  Less 
than  20  ft  of  bank  erosion  was  reported  at  the  upstream  end  of  the  reach  near 
the  bridges  and  less  than  50  ft  at  the  downstream  end  near  Ward  Canyon.  The 
existing  masonry  floodwall  at  the  bridges  sustained  some  damage  during  the 
1983  floods.  Repair  work  observed  during  a  May  1988  field  reconnaissance 
indicated  that  damage  was  caused  by  floodflows  returning  to  the  channel  rather 
than  by  undercutting  of  the  toe. 
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THALWEG  ELEVATION 


THALWEG  ELEVATION  IN  FEET  THALWEG  ELEVATION  IN  FEET 


Figure  6.  Sensitivity  of  calculated  thalweg  to  sediment 
inflow  at  end  of  standard  project  flood 


Figure  7.  Sensitivity  of  calculated  thalweg  to  bed  material 
gradation  at  peak  standard  project  flood 
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30.  Riffle  and  pool  sequences,  characteristic  of  gravel-bed  rivers, 
typically  have  regular  spacings  with  impinging  flows  at  pools  on  alternating 
sides  of  the  channel.  The  San  Francisco  River,  through  Clifton,  even  with  its 
width  constraints,  has  developed  a  typical  riffle-pool  sequence  (Plate  17). 
Aerial  photos  taken  after  the  1983  flood  reveal  a  consistent  pattern  upstream 
of  the  railroad  bridge  with  pools  alternating  from  one  side  of  the  channel  to 
the  other,  with  spacings  of  5  to  8  times  the  channel  width  or  about  1,000  ft. 
Downstream  from  Clifton,  where  the  channel  is  no  longer  encroached  on  by  de¬ 
velopment,  the  natural  channel  also  has  a  consistent  sequence  of  alternating 
pools  and  riffles,  with  spacings  of  about  five  channel  widths  or  1,800  ft. 

The  consistent  alternating  sequence  upstream  and  downstream  of  the  proposed 
leveed  reach  is  discontinuous  in  the  leveed  reach  itself.  If  the  upstream 
sequence  were  to  extend  downstream  past  the  highway  bridge,  the  flow  would 
impinge  against  the  existing  left-descending  bank  near  sta  552+00.  This  is 
the  area  where  the  continuity  of  the  slag  deposit  is  uncertain,  and  the  pro¬ 
posed  levee  alignment  sets  on  top  of  the  inferred  slag  deposit  alignment. 
However,  the  existing  thread  of  the  channel  does  not  follow  the  upstream 
pattern  and  does  not  impinge  on  the  left-descending  bank  until  further  down¬ 
stream  where  the  slag  deposit  is  exposed,  and  the  levee  alignment  is  set  back 
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100  ft.  Comparison  of  the  1939  topographic  map  and  the  1983  aerial  photos 
indicates  that  the  low-flow  channel  has  moved  away  from  the  left  descending 
bank  adjacent  to  the  proposed  levee  between  river  stations  544+00  and  552+00. 
This  favorable  condition  is  attributed  to  longitudinal  progression  of  the 
gravel  bar  along  the  left  descending  side  of  the  channel.  Current  plans  to 
excavate  this  gravel  bar  could  result  in  unfavorable  realignment  of  the  low- 
flow  channel  toward  the  left  bank. 

31.  Historical  survey  data  were  reviewed  to  evaluate  riverbed  stability 
over  time.  Historical  profile  and  cross-section  data  were  published  in  the 
Planning  Assistance  Study  (USAED,  Los  Angeles,  1979)  and  are  shown  in 
Plates  18-22.  About  4  ft  of  degradation  occurred  in  the  low-flow  channel  at 
the  downstream  end  of  the  project  reach  between  1939  and  1977.  This  condition 
was  confirmed  during  site  reconnaissance  in  May  1988.  The  degradation  has 
undercut  the  adjacent  slag  deposit  and  exposed  a  utility  pipe  crossing  the 
river.  Upstream  migration  of  the  headcut  has  apparently  been  arrested  by  a 
collection  of  large  rocks  near  sta  539+00.  These  boulders,  some  reaching 
20  ft  in  diameter,  extend  from  the  left-descending  bank  across  the  low-flow 
channel,  and  about  100  ft  onto  the  adjacent  gravel  bar.  There  is  no  observ¬ 
able  bedrock  at  this  location;  the  boulders  are  reported  to  be  remnants  of  a 
rock  promontory  from  the  right  bank  that  was  blasted  away  in  1979  (USAED, 

Los  Angeles,  1988).  Structural  measures  such  as  this  may  be  required  to  pre¬ 
vent  future  degradation  of  the  low-flow  channel  and  ensure  a  thalweg  elevation 
above  the  slag  deposit  toe. 


Regression  Analysis 

32.  A  river  that  erodes  its  boundary  will  form  a  channel  that  effec¬ 
tively  carries  its  water  and  sediment  load.  Theories  of  equilibrium  channel 
formation  are  based  on  the  premise  that  hydraulic  parameters  can  be  deter¬ 
mined,  on  the  average,  as  a  function  of  an  independent  variable,  usually  dis¬ 
charge.  Regression  equations  have  been  developed  by  a  number  of  investigators 
for  width,  depth,  velocity,  resistance,  and  slope.  Of  these,  predictors  for 
width  are  the  most  reliable. 

33.  Hydraulic  geometry  regression  equations  were  determined  by  Leopold 
and  Maddock  (1953)  for  selected  rivers  in  the  midwestern  United  States.  The 
generalized  equation  for  width  is  expressed: 
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where : 


w  =  width 

Q  =  characteristic  discharge 

a  =  coefficient  characteristic  of  the  particular  river 
b  =  exponent  determined  to  generally  be  a  constant 
Potential  for  bank  erosion  adjacent  to  the  proposed  levee  at  Clifton  was  eval¬ 
uated  by  calculating  predicted  widths  using  the  regression  equations  of  North¬ 
west  Hydraulic  Consultants  (NHC)  and  Williams  (1978). 

34.  NHC  (1984)  proposed  that  hydraulic  geometry  could  be  used  for  pre¬ 
liminary  design  of  stable  channels.  They  came  up  with  an  exponent  b  of  0.50 
for  gravel  bed  streams.  NHC  also  determined  that  the  coefficient  a  could  be 
determined  depending  on  the  nature  of  the  bank  material.  Banks  with  high  re¬ 
sistance  (riprap,  cobbles  and  boulders,  or  stiff  clay)  are  assigned  a  coeffi¬ 
cient  of  1.6,  moderate  resistance  (coarse  gravel,  medium  clay)  are  assigned  a 
coefficient  of  2.1,  and  low  resistance  (sand  and  silt,  soft  clay)  are  assigned 
a  coefficient  of  2.7.  In  the  NHC  procedure,  it  is  recommended  that  either  the 
10-yr-frequency  peak  discharge  or  50  percent  of  the  design  discharge,  which¬ 
ever  is  larger,  be  used  for  Q  in  the  regression  equation.  Using  the  NHC 
procedure  with  a  coefficient  of  2.1,  widths  of  565  and  606  ft  were  calculated 
for  discharges  of  72,500  and  83,500  cfs,  respectively  (50%  of  design  and  SPF 
discharges) . 

35.  At-a-station  hydraulic  geometry  regression  equations  consider  tem¬ 
poral  variations  in  channel  geometry  as  discharge  fluctuates  at  a  cross  sec¬ 
tion.  The  exponent  b  in  equations  for  width  are  significantly  affected  by 
bank  cohesion.  A  river  with  near-vertical,  cohesive  banks  would  have  an  ex¬ 
ponent  near  zero  because  the  width  would  essentially  be  constant  with  changes 
in  discharge.  Based  on  empirical  data  from  various  types  of  channels, 

Williams  (1978)  came  up  with  an  exponent  of  0.40  for  sandy-gravelly  rivers 
with  one  cohesive  bank  and  one  noncohesive  bank.  According  to  Williams,  this 
result  compares  to  a  theoretical  value  of  0.45  using  Langbein's  (1965)  minimum 
variance  theory.  Coefficients  for  Williams'  equation  were  determined  for  the 
San  Francisco  River  using  an  aerial  photo  dated  7  October  1983.  An  average 
channel  width  of  110  ft  in  the  study  reach  was  measured  from  the  photos. 

Using  the  mean  daily  flow  (1110  cfs)  for  7  October  1983,  a  coefficient  of  6.7 
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was  calculated.  Using  the  channel  trace  cleared  of  vegetation  as  Indicative 
of  channel  width  at  the  flood  peak,  an  average  channel  width  of  340  ft  was 
determined  from  the  aerial  photos  for  the  flood  that  occurred  on  2  October 
1983.  Using  the  mean  daily  flow  (52,200  cfs)  for  the  discharge  in  the  hydrau¬ 
lic  geometry  equation,  a  value  of  4.4  was  calculated  for  coefficient  a  . 

Using  the  peak  flow  of  90,900  cfs,  a  value  of  3.5  was  calculated  for  a  .  The 
hydraulic  geometry  equation  using  an  exponent  of  0.40  and  a  coefficient  of  3.5 
yielded  a  predicted  width  of  406  ft  and  430  ft  for  peak  discharges  of 
145,000  cfs  and  167,000  cfs,  respectively.  With  a  coefficient  of  6.7,  pre¬ 
dicted  widths  of  777  ft  and  873  ft  were  calculated  for  peak  discharges  of 
145,000  cfs  and  167,000  cfs,  respectively. 

36.  Proposed  widths  between  the  existing  rock  right  bank  and  the  pro¬ 
posed  levee  varies  between  300  and  500  ft.  This  analysis  indicates  that,  un¬ 
restrained,  there  would  be  a  tendency  for  the  river  to  widen  past  the  proposed 
levee  alignment.  This  regression  analysis  does  not  account  for  the  erosion- 
resistant  slag  deposit  adjacent  to  the  left  descending  bank.  In  addition,  the 
application  of  these  generalized  equations  to  flashy  streams  like  the  San 
Francisco  River  is  questionable.  It  can  be  concluded  however  that  the  integ¬ 
rity  of  the  slag  deposit  is  important  because  it  arrests  the  river's  tendency 
to  widen  at  high  flows. 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


37.  Results  of  the  numerical  model  study  indicated  that  the  channel  bed 
elevation  would  vary  significantly  during  the  passage  of  a  major  flood.  Sig¬ 
nificant  scour  would  occur  downstream  from  the  Highway  666  bridge  where  exist¬ 
ing  bank  protection  is  provided  by  a  masonry  wall  with  unknown  toe  protection, 
and  the  proposed  levee  toe  would  have  less  than  the  100-ft  setback.  Degrada¬ 
tion  in  the  lower  portion  of  the  project  reach,  during  the  passage  of  a  major 
flood,  is  prevented  by  a  backwater  condition  created  by  a  downstream  constric¬ 
tion.  However,  degradation  in  this  lower  reach  occurs  during  low  flows  when 
the  backwater  condition  does  not  exist. 

38.  The  review  of  historical  data,  including  aerial  photos,  surveyed 
profiles,  and  eye  witness  accounts,  indicates  that  there  is  some  potential  for 
bank  erosion  in  the  reach  downstream  from  the  bridge  at  sta  566+80  during  high 
flows  and  in  the  lower  reaches  during  lower  flows. 

39.  Hydraulic  geometry  regression  analysis  indicates  there  would  be  a 
tendency  for  channel  widening  during  flood  flows. 

40.  The  existing  masonry  wall,  downstream  from  the  bridges,  is  in  a 

zone  of  high  bank  erosion  and  degradation  potential.  It  is  recommended  that 

the  existing  toe  elevation  be  identified,  and  if  necessary,  additional  toe 

* 

protection  provided.  Numerical  model  study  results  indicate  scour  depths  on 
the  order  of  5  ft  with  existing  channel  geometry.  However,  these  calculations 
are  qualitative.  Engineering  judgment  must  be  relied  on  heavily  to  determine 
an  adequate  depth  for  toe  protection.  Downstream  from  the  masonry  wall,  the 
levee  alignment  lies  cn  top  of  the  inferred  slag  deposit.  In  this  reach,  the 
slag  must  be  considered  an  integral  part  of  the  levee  design  and  toe  protec¬ 
tion  is  recommended  to  insure  the  integrity  of  the  slag  deposit.  Historical 
data  indicate  that  the  gravel  bar  in  this  location  has  existed  for  a  long 
time,  and  provides  some  toe  protection  for  bank  improvements.  However,  if 
removed,  the  flow  could  redirect  itself  and  move  toward  the  left-descending 
bank.  It  is  recommended  that  the  gravel  bar  be  allowed  to  remain  in  place. 

41.  Downstream  from  sta  550+00,  the  existing  slag  deposit  along  the 
left  descending  bank  is  critical  to  the  integrity  of  the  proposed  levee.  In 
areas  where  the  slag  deposit  is  being  undercut,  remedial  efforts  are  recom¬ 
mended  to  halt  undermining  of  the  existing  toe.  This  could  be  accomplished  by 
dumped  stone  along  the  toe  or  by  grade  control  structures.  The  longitudinal 
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continuity  of  the  slag  deposit  is  essential  to  the  integrity  of  the  levee. 
Its  existence  should  be  confirmed. 


42.  Downstream  from  sta  550+00,  setback  distance  for  the  levee  is  not 
especially  critical  with  respect  to  bank  erosion.  As  long  as  the  slag  deposit 
remains  in  place,  the  levee  toe  will  be  protected.  Some  slumping  or  breaking 
of  the  slag  deposit  during  a  single  flood  event  would  not  affect  the  proposed 
levee  toe.  The  mass  of  the  deposit  would  be  sufficient  to  continue  to  provide 
toe  protection  for  the  duration  of  the  flood.  However,  it  is  recommended  that 
the  integrity  of  the  slag  deposit  be  included  as  a  maintenance  requirement  for 
the  project.  The  levee  should  be  set  back  far  enough  so  that  the  levee  toe  is 
not  structurally  effected  by  the  presence  of  the  slag  deposit.  Setback  dis¬ 
tance  can  be  determined  based  on  real  estate,  structural,  geotechnical,  and 
construction  considerations. 
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